Dynamical mean field theory for manganites 
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Doped and undoped manganites are modeled by the coupling between itinerant Cg electrons and 
static t2g spins, the Jahn- Teller and breathing phonon modes, and the Coulomb interaction. We pro- 
vide for a careful estimate of all parameters and solve the corresponding Hamiltonian by dynamical 
mean field theory. Our results for the one-electron spectrum, the optical conductivity, the dynamic 
and static lattice distortion, as well as the Curie temperature show the importance of all of the above 
ingredients for a realistic calculation as well as for describing the unusual dynamical properties of 
manganites including the insulating parent compound and the insulating-like paramagnetic state of 
doped manganites. 

PACS numbers: 71.27.-|-a, 71.10.Fd, 75.47.Gk 
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I. INTRODUCTION 

Manganites have attracted intensive interest dur- 
ing the last decades due to their extraordinary 
properties including the colossal magnetoresistance 
(CMR)ii^i^. These materials have the chemical composi- 
tion Ti-jjDajMnOa, where T is a trivalent rare earth ion 
(T = La, Pr, Nd, . . .) and D is a divalent alkali ion (D 
— Ca, Sr, . . .) and crystallize in a cubic perovskite struc- 
ture with a possible distortion at low temperatures, see 
Fig.[Tl which albeit can be distorted at low temperatures. 
Soon after CMR was discovered, various phase diagrams 
as a function of temperature and magnetic field or dop- 
ing concentration x were established'^'^'^i^'^'^ i^^'^-'^i"'^^ ; also 
see the review articles Refs. [T3lll^lT5lll6l . For the parent 
compound LaMnOa, an insulator-to- metal transition was 
found upon applying pressureii; and a chargc/orbital- 
ordered phase has been reported in a large number of per- 
ovskite (Ti_a;D^Mn03) and layered (Ti_2;Di+::cMn04) 
manganites, depending on the effective bandwidth and 
the quenched disorde r^^'^^ . The dynamical properties 
of the paramagnetic insulating are very unusual as is 
reflected in a spectral function A{lu) with a very low 
spectral weight at the Fermi level Ep irrespectively of 
X , as indicated by photoemission and X-ray absorption 
experiments2£ii2ii22ii^. Similarly, the optical conductiv- 
ity 17(0;) shows a very low spectral weight up to an en- 
ergy scale of ~ 1 eV— Besides, the ferromagnetic 
metallic phase is an atypical (bad) metal^i. 

A physical understanding of these properties is diffi- 
cult due to the internal complexity resulting from the 
interplay between charge, spin, orbital and lattice de- 
grees of freedo m^^i^^ . A basis ingredient of a theoreti- 
cal description is the separation of the electrons within 
the five d orbitals into a localized t2g spin of length 
|S| — 3/2 and n — I ~ x itinerant eg electrons, cou- 
pled to the t2g spin by Hund's exchange. This the ba- 
sis of the so-called "double exchange" ^^i"^^ which led to 
the (ferromagnetic) Kondo lattice model^^ and explains 
ferromagnetism in doped manganites and the charge- 
ordered phase at x=0.5^. A spin-canted state was sug- 
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FIG. 1: Sketch of the cubic perovskite unit cell for mangan- 
ites. The dashed lines indicate the basic MnOe octahedron. 



gested later'^-. However, this double exchange model- 
ing disagrees with the experiment in many aspects in- 
cluding the CMR. Hence, the importance of the Jahn- 
Teller phonon modes and their couling to the eg electrons 
was stressed"^^ and studied^'^'^'^. However, while it 
can describe insulating-like behavior for large electron- 
phonon coupling when electrons are trapped as lattice 
polarons^l, the Kondo lattice model extended by Jahn- 
Teller phonons still fails to produce a large magnetore- 
sistance at finite doping. Another important ingredient 
for the physics of manganites is the local Coulomb in- 
teraction between the eg electrons as was pointed out 
in""^ and in realisitic ab initio calculations^i^i^i^. Of 
these ab initio calculations, the combination of the local 
density approximation (LDA) and dynamical mean field 
theory {DMFT)^d^dL^^dl>^ has the least tendency to 
overestimate the formation of an insulating state since 
the DMFT electron dynamics also avoids the cost of dou- 
ble occupations in the paramagnetic metallic state. The 
corresponding LDA-f DMT results for manganites^^ show 
that only the combined localization effect of Coulomb in- 
teraction and (static) Jahn- Teller distortion makes the 
undoped parent compound an insulator, and leads to a 
correct desription of the pressure induced insulator-to- 
metal transition^. Similarly, for the doped compound, 
both mechanisms work together in localizing the eg elec- 
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trons giving rise to a proper description of the unusual 
dynamical properties of doped manganites and the large 
CMR in a wide range of doping^. 

In this paper, we present a detailed theoretical in- 
vestigation of manganites using DMFT— i^i^ to solve 
a realistic model of manganites which supplements the 
Kondo-lattice model by Jahn- Teller and breathing mode 
phonons as well as by the local Coulomb interaction be- 
tween the eg electrons; first results have been presented 
in^i also see the recent DMFT calculation The pa- 
per is organized as following: Section |TT] introduces the 
Hamiltonian and discusses how realistic parameters are 
choosen. Section IIIII is devoted to the undoped parent 
compound LaMnOs focusing in Section UlI Al on how the 
Jahn- Teller phonons, for a cubic lattice, give rise to a 
very similar spectrum as that of LDA-I-DMFT for the 
experimental static Jahn- Teller distortion. This is sup- 
plemented by Section FlII BI studving the symmetry break- 
ing towards a static Jahn- Teller distortion. Section IIVI 
presents results for doped mangaites including a detailed 
study of the paramgnetic insulating-like phase in Sec- 
tion IIV Al and of the phase transition towards ferromag- 
netism in Section flV Bl Section IVl includes the effect of 
the breathing mode with an estimate of the Coulomb 
repulsion and the electron-phonon coupling strength in 
Section IV Al and the corresponding optical conductivity 
in Section IV Bl Finally, Section IVII gives a summary of 
the results obtained. 



II. REALISTIC MODEL AND PARAMETERS 

For the realistic microscopic modeling of manganites, 
we employ the Hamiltonian 



d,j.2_y2 and d^z'^^.r2 indexed by /i (and v) 
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Here, cj^^ and Ci^a are the fermionic creation and anni- 
hilation operators for electrons at site i within Cg-orbital 
fi and with spin cr; s^^ is the corresponding spin operator 



with Pauli matrices r. 



The first line of Hamiltonian ^ forms the ferromag- 
netic Kondo lattice model with a coupling J' of the Cg 
spin to the localized t2g spin S^. For the cubic lattice, the 
hopping elements can be, to a good approximation, 
restricted to neighboring sites i and j in the x, y and z 
direction with an orbital matrix for the two (cg) orbitals 
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The constant to = W/6 = 0.6 eV was estimated from 
LDA which gives a bandwidth W = 3.6 eV for cubic 
LaMnOg^. From the 2.7 eV splitting between spin-up 
and spin-down Cg bands in the corresponding ferromag- 
netic LDA calculatation on the other hand, the value of 
coupling to the t2g spins was estimated 21/18^1 = 2.7eV. 

The second line of Hamiltonian ([1]) describes the 
Coulomb repusion of two Cg electrons on the same lat- 
tice site i, consisting of the intra-orbital Coulomb repul- 
sion U, the inter-orbital repulsion U' and the Hund's rule 
energy gain J for two parallel eg spins. The Coulomb in- 
teraction can be estimated from the experimental spec- 
trum, combining photoemission spectroscopy (PES) and 
x-ray absorption spectroscopy (XAS)^. Accounting for 
the crystal field splitting between eg states and t2g states 
of 1 — 2 eV, we obtain an average Coulomb interaction 
U = 3 — 4: cV, which is also supported by spectral ellip- 
sometry for LaMnOg^. Since U' = J7 — 2J by symmetry, 
we have the relation 

U^^^^-''^''-''^U-'-J. (4) 
3 3 ^ 

Taking a value J = 0.75 eV, which is slightly smaller than 
the value for the five band model obtained by constrained 
LDAS and which agrees with our own estimate for J7, 
all Coulomb interactions are determined. Unless noted 
otherwise, we employ U' = 3.5 eV. 

The third line is the Hamiltonian for the three most 
important phonon modes Qa, i.e., the breathing mode Qi 
and the two Jahn- Teller modes Q2,3 as illustrated in Fig. 
[21 The only free parameter of this part of the Hamilto- 
nian II]) are the three phonon frequencies which, follow- 
ing Ref. [ssl . can be obtained from the Raman spectrum 
for LaMnOa. The Raman spectrum gives very similar 
values for the Jahn- Teller mode frequency Ojt = 0.07 eV 
and that of the breathing mode flbr = 0.08 eV so that we 
take a unique value Q — 0.07 eV in this paper. 

Finally, the fourth line is the electron-phonon coupling 
with the breathing mode coupling to the electron density 
and the two Jahn- Teller modes coupling to the difference 
in eg occuption. Avoiding to use too many parameters, 
we take, as for the frequency, a unique coupling strength 
g. For the Jahn- Teller modes it is related to the static 
Jahn- Teller energy Ejt — 5^/257^, the ground state en- 
ergy if only the Jahn- Teller coupling and phonon energy 
is present. One can try to determine the Jahn- Teller cou- 
pling stength g from the lattice distortion. A distortion of 
0.1 A which is consistent with some LDA calculations^!^ 
leads to Ejt = 0.25 eV and hence g = 0.05 eV^/^ for 
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FIG. 2; Vibration of the oxygen octabedra around the man- 
ganese ion for the breathing mode Qi and the two Jahn-Teller 
modes Q2 and Q3. 



il ~ 0.07 eV. However, recent x-ray powder diffraction 
and neutron powder diffraction measurements found a 
much larger distortion^25.^ which would result in unrealis- 
tically large values for g or frequencies at odds with the 
Raman frequencies. We hence attribute these larger dis- 
tortions to cooperative lattice effects and the quadratic 
vibronic coupling to the electronic degrees of freedom. 
Possibly also enhancement effects due to electronic cor- 
relations play a role. As a consequency the precise value 
of g is an open issue and we have hence done calculations 
for various values oi g. In section IV Al we provide for a 
new estimate of the two parameters with the biggest un- 
certainty, i.e., U and g on the basis of ther experimental 
gap of undoped LaMnOa and the resistitvity for doped 
LaMnOa. 

For the folowing results, Hamiltonian ([T]) is solved 
using DMFT^i^i^ with Hirsch-Fye^^ Quantum 
Monte Carlo (QMC) simulations supplemented 
by the Blankenbecler-Scalapino-Sugar algorithm^^ 
for Holstein phonons. To this end, the inverse 
temperature (3 is discretized into L time slices 
Ti = {l- 1) At, = 1, • • • , L 1) of size At = 0.25 eV'^ 
and the t2g spin is assumed to be classical. Let us briefly 
discuss some aspects of the phonon fields since these are 
less commonly simulated in DMFT. The phonon field 
can be described by a classical field with boundary 
condition (pL+i = </>!—■ In the effective action, the 
kinetic and potential energy of the phonons [third line 
in ll])] hence become 
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We employ local updates for one time slice / and global 
updates for all time slices, both having the form (j)^ = 
(j^i -\- {x — l/2)5(j) with a random number x £ [0,1) 
and a properly choosen amplitude 5(j) different for lo- 
cal and global updats. For low temperatures T and 
large electron-phonon couplings g, the probability to go 
through the QMC sampling from a region with large lat- 
tice distortion (occupied with one eg electron) to a re- 
gion with small lattice distortion without Cg electron be- 
comes very low. This is quite similar to the Hubbard 
model where for large U and small T, the tranisition 
between configurations with predominately spin-up and 
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FIG. 3: Local spectral density in the paramagnetic phase 
of LaMnOs calculated for g = 0.05 eV^^^ U' = 3.5 eV and 
/3 = 30 eV~^. All other parameters are fixed as discussed in 
Section]^ The breathing mode phonon will only be included 
later in Sec. [3 



spin-down becomes very rare. Fortunately, the physics 
in both situations is rather simple and a proper aver- 
aging between spin-up and spin-down respectively small 
and large lattice distortion is easily possible. In the latter 
(manganites) case, the probability V = x for unoccupied, 
undistorted sites is physical meaningful and also obtained 
by an extrapolation from higher temperatures. 

From the DMFT one-electron spectral function the op- 
tical conductivity is calculated from the simple bubble 
diagram given by two Green functions multiplied by two 
group velocities (derivates of the dispersion relation) . For 
the magnetic transition temperature Tc, the two particle 
Green function and from that the inverse susceptibiliy 
crossing zero at were calculated. 



III. INSULATING PARENT COMPOUND 

A. Paramagnetic phase 

We start our study with the undoped parent com- 
pound LaMnOs, which for a static Jahn-Teller distor- 
tion was studied in Ref4^ using LDA-hDMFT. Fig. [3] 
shows the spin- and orbital-averaged spectral density 
A{lj) ~ S^cr "^^'^'^('^) paramagnetic phase 

with the Green functions G^'^(uj) obtained for real fre- 
quencies uj by the maximum entropy method^. It shows 
a three preak structure very similar to the previous 
LDA-fDMFT resultsi^ even though no static Jahn-Teller 
distortion is present. The reason for this is that the lat- 
tice dynamics leads to a similar (averaged) lattice dis- 
tortion as in the static case, see Fig. 2] At a given time 
shot (or for a given QMC configuration), the two Cg lev- 
els are Jahn-Teller split with the precise decomposition 
of the lower (and higher) lying orbital depending on the 
sign and admixture of the Q2 and Q3 Jahn-Teller modes. 
The first peak around -1.3 eV corresponds to electronic 




Lattice distortion (A) 

FIG. 4: Probability distribution (in arbitrary unit) of the 
lattice distortion for the Jahn-Teller coupling g = 0.05 eV'^^"^ , 
Coulomb repulsion U' = 3.5 eV and inverse temperature (5 = 
16 eV-^ 
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FIG. 5: Local spectral density for the alternating orbital 
ordering calculated by DMFT (QMC). The parameters are 
U' = 3.5 eV, T = 0.05 eV and g = 0.05 eV^/^ 



states occupying the lower of the two levels, and the 
second one around +1.3eV to adding an electron in the 
higher lying orbital. The third peak around +4.5 eV 
corrsponds to excitations to states with antiparallel t2g 
spin, whereas the first two peaks have parallel t2g spin. In 
Fig. [31 there is still some spectral weight in the gap which 
may stem from the phonon sideband^**. The energy gap 
can hence only be obtained approximately. It is compa- 
rable to the experimental gap of about 1 eV— i^i. The 
consistency between the previous LDA+DMFT and the 
present DMFT model calculations supports our realistic 
microscopic model ([l} for describing the electronic be- 
havior of LaMnOa and demonstrates once again that the 
insulating ground state in LaMnOa at ambient conditions 
results from the combination of the Coulomb interaction 
and the Jahn-Teller coupling in addition to the Hund's 
coupling between the and t2g spins. 

Fig. |4] plots the probability distribution P{Q) for a 
given lattice distortion in the range [QtQ + dQ]. Note 
that the lattice distortion are converted into units of A 
by multiplying a factor of h/^/M where M is the mass 
of oxygen atom and that for a finite At the two modes 
are not exactly symmetric. 

We see that P{Q) has one broad peak located at 
around 
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0.167 A, 



(6) 



which corresponds roug hly to Q* = g/^"^ w 0.15 A ob- 
tained for a single-site model. As we have discussed, the 
discrepancy from the neutron experiment^*' is probably 
due to the higher order quadratic vibronic coupling and 
the cooperative effect between the adjacent MnOg Oc- 
tahedra that share a common oxygen atom^. A recent 
LDA-I-U calculation has taken into account these effects 
and produced the correct experimental results^. Clus- 
ter extensions are required for further realistic DMFT 
calculations. 



B. Structural transition 

Experimentally a structural first-order phase transi- 
tion is observd at Tqo ~ 740 K with an abrupt vol- 
ume contraction^. The lattice is nearly cubic above 
Too but has a strongly distorted orthorhombic structure 
due to the static Jahn-Teller distortion below Too- The 
structural transition is accompanied by an orbital order- 
disorder transition. The low temperature phase shows a 
staggered ordering of the d3x2_^2 and d^y2_^2 orbitals in 
the a-b plane which repeats itself along the c-direction. 

With a slight modification of the single-impurity 
DMFT (QMC) algorithm, we can study a hypotheti- 
cal antiferromagnetic orbital ordering on an AB lattice 
and draw some conclusions about the structural transi- 
tion in LaMnOs. To study the alternating orbital or- 
dering, the cubic lattice is separated into two sublattices 
A and B with oposite behavior of the two orbitals: 
S^2)('^) = -Due to the symmetry, we only need 

to take care of a single lattice site of either type in DMFT. 

Fig. [5] shows the local spectral densities of two Cg or- 
bitals for U' = 3.5 eV, T = 0.05 eV and g = 0.05 eV^/^. 
At this temperature, the orbital symmetry is strongly 
broken: One orbital is occupied with almost one electron, 
while the other orbital is only slightly occupied. Fig. [S] 
also indicates the contributions from different spin and 
orbital components. A direct consequence of the orbital 
order is, due to the Jahn-Teller coupling, a correspond- 
ing lattice distortion. Averaged over the two orbital, the 
spectrum is actually very similar to that of Fig. [3] with 
the two peaks at — 1 eV and 2 eV stemming from the spin 
states parallel to the local t2g spins. These two peaks 
are split by a combination of Jahn-Teller coupling and 
Coulomb interaction. 

Fig. [B] plots the orbital polarization P = \ ^„[nf^ — 
n^^)|, i.e., the difference of orbital occupations rt^g.; ^ 
a function of bandwidth and temperature. For W = 
3.6 eV, we find a finite orbital polarization below T w 
725 K, in agreement with the experimental result of the 
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FIG. 6: Orbital polarization vs. temperature at fixed W — 
3.6 eV (squares) and vs. bandwidth at fixed T = 0.05 eV (cir- 
cles). The Jahn-Teller coupling is g = 0.05 eV"^^^ and the 
inter-orbital Coulomb repulsion U' — 3.5 eV. 



structural transition temperature Tqo = 740 K^. If we 
fix the temperature T = 0.05 eV and var the bandwidth 
as under pressure, the orbital polarization is almost a 
constant for the bandwidth W < 4.8 eV and then de- 
creases with increasing bandwidth until it is reduced to 
zero a.t W = 6.0 eV. This behavior reflects the nature 
of the bandwidth-control metal-insulator transition in 
LaMnOa. Here, W = 4.8 eV corresponds to the criti- 
cal bandwidth where the split minority and majority Cg 
bands start to overlap; and LaMnOa becomes metallic. 
Hence, for a larger bandwidth, the orbital polarization 
decreases rapidly with increasing bandwidth. 

The bandwidth W = 6.0 eV marks a second transition 
where the orbital polarization and the Jahn-Teller distor- 
tion are completely suppressed. Although different from 
Ref. [13, we can still identify three distinct regimes at low 
temperature: (i) an insulating phase with orbital order- 
ing and static Jahn-Teller distortion below W « 4.8 eV 
(or P/j\/ — 32GPa); (ii) an intermediate metallic phase 
with orbital ordering and static Jahn-Teller distortion be- 
low W ~ 6 eV; (iii) a metallic phase with orbital symme- 
try and dynamic Jahn-Teller distortion above W ^ 6 eV. 



IV. DOPED MANAGANITES 
A. Paramagnetic insulating state 

Let us now turn to the doped manganites with their 
extraorinary properites such as the colossal magnetore- 
sistance. For simplicity, we first neglect the breathing 
phonon mode. We start by plotting the probability dis- 
tribution of the lattice distortion for n — 0.8 {x — 0.2) 
at different temperatures and couplings in Fig. \7\ The 
Coulomb interaction is fixed to U' = 3.5 eV throught 
Section IIVI In contrast to the single peak distribution 
obtained for the undoped case in Fig. 21 we find two 
peaks in the distribution function at large Jahn-Teller 
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FIG. 7: Probability distribution of the lattice distortion at 
n = 0.8 Cg electrons per site (i.e., doping x — 0.2) for (a) 
g = O.lOeV^/2 and (3 = 8 and 16 eV"^ and (b) (3 = 16 eV"^ 
ajid g = 0.05, 0.10, and 0.15 eV^^^ 



couplings for doped systems. The large peak located 
at about g/Vl^ corresponds to the large lattice distor- 
tion due to the Jahn-Teller coupling (and corresponds to 
the single peak without doping), while the small peak 
stems from the quantum and thermal fluctuation of the 
MnOg octahedra which due to the doping are not occu- 
pied with an electron (and hence not strongly split by the 
electron-phonon coupling g). While at a large coupling 
strength g ^ 0.1, the two peaks are well separated, they 
merge into a single structure at weak coupling such as 
for g ^ 0.05 eV^/^ in Fig. [H 

The connection to the electronic spectrum is shown in 
Fig. [51 For clarity, we separate the spectrum into two 
parts originated from the large and small lattice distor- 
tions, respectively. For the large lattice distortion, the 
corresponding spectral density locates well below or far 
above the Fermi energy. The low energy part can be iden- 
tified as localized eg electrons or more precisely the tran- 
sition from this state to states without Cg electron. The 
high energy part stems from the states which are pushed 
up by the Hund's rule coupling, the Jahn-Teller splitting 
and the Coulomb repulsion. For the small lattice distor- 
tion, most of the spectral density locates slightly above 
the Fermi energy. They stem from the " undistorted" un- 
occupied states which are also called "midgap states" in 
the literature^ML. 

The formation of lattice polarons provides for the ba- 
sic physics of doped manganites. The electron spec- 
tral density can be seen as a combination of polaron 
states well below the Fermi energy and the midgap states 
above the Fermi energy. As shown in Fig. [51 this results 
in the strong suppression of the spectral weight at the 
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FIG. 8: Probability distribution of the lattice distortion and 
the corresponding electron spectral densities for the Jahn- 
Teller coupling g = 0.10 eV^/^ (left panel) and 0.15 (right 
panel) at inverse temperature (3 = 16eV~^ and n = 0.8 eg 
electrons per site. The electron spectral densities are sep- 
arated into two parts with large and small lattice distor- 
tions. The arrows indicate the separation between peaks. The 
phonon distribution at g = 0.10 eV'^''^ without the Coulomb 
interaction is also shown for comparison. 



Fermi energy and gives rise to a large energy gap for 
g = 0.15 eV'^/^ and a pseudo-gap for smaller Jahn- Teller 
couplings. 

Let us discuss the role of the on-site Coulomb interac- 
tion in the formation of the pseudo gap and hence the 
insulating-like paramagnetic state of manganites. Fig. [S] 
compares the spectral densities for g — 0.10 eV^^^ with 
and without the Coulomb interaction; Fig. [S] the cor- 
responding probability distribution of the lattice distor- 
tion. For U = 5eV {U' = 3.5 eV), some spectral weight 
is pushed away from the Fermi energy to higher ener- 
gies and the pseudo-gap is strongly enhanced by the 
Coulomb interaction. Also the lattice distortion shows 
a more pronounced separation of the two peaks which 
indicates the enhancement of the polaron formation by 
the Coulomb interaction. This confirms the important 
role of the Coulomb interaction. 

Let us now consider the gradual swiching on of the 
Jahn- Teller couplings g at fixed /? = 16eV~^ and n = 
0.8 shown in Fig. [TOl Without Jahn- Teller coupling, the 
spectrum has a broad quasiparticle peak at the Fermi 
energy and the system is metallic. With increasing g, 
the quasiparticle peak is gradually suppressed and the 
spectral weight at the Fermi energy is reduced, leaving 
a dip in the spectrum for intermediate coupling and an 
energy gap for strong coupling. The low (high) energy 
Hubbard band also shifts towards lower (higher) energies. 
An additional peak shows up at w = 1 — 2 eV, ascribed to 
the midgap states with small dynamic lattice distortion. 

Turning to the doping dependence at fixed P = 
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FIG. 9: Paramagnetic electronic spectral density at n =0.8, 
f3 = 16 eV"^ and g = 0.10 eV^/^ comparing U = OeV and 
U = 5eV (U' = 3.5 eVj. 
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FIG. 10; Electron spectral density in the paramagnetic phase 
for different couplings g = 0, 0.05, 0.10, and 0.15 eV^/^. The 
parameters are (3 = lGeV~^ and n = 0.8. The quasiparticle 
peak is suppressed at large g, giving the pseudo-gap behavior 
for intermediate coupling and a large energy gap for strong 
coupling. 

16 eV"^ and g = 0.10 eV^/^ ^ig. [Ill we see that the 
system is a good insulator with a large energy gap of 
about 2.3 eV without doping. This is larger than the ex- 
perimental gap of 1 eV and hence suggests that a smaller 
Jahn- Teller coupling such as g — 0.05 eV'^/^ used in Sec- 
tion mil is more realistic. But for the qualitative dis- 
cussion here, we continue with g = 0.1 eV^^^ since for 
this larger Jahn- Teller coupling the different features of 
the spectrum can be better identified. A small doping 
a; = 0.1 introduces some midgap states just above the 
Fermi energy and changes the large energy gap into a 
pseudo-gap at the Fermi energy. If the doping is large 
enough, the midgap states can dominate at the Fermi 
energy so that the pseudo-gap is completely filled at 
n = 0.3. However, due to the strong phonon and spin 
scattering, the quasiparticle peak is still damped (broad- 
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FIG. 11: Paramagnetic spectral density for different electron 
occupations n —1.0, 0.9, 0.8, 0.5, and 0.3 at (3 = 16eV~^ and 
g = 0.10 eV^/^ 
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doping, X 

FIG. 13: Curie temperature as a function of doping x — l~n 
for g = 0.10 eV^'''^ 
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FIG. 12: Optical conductivity in the paramagnetic phase at 
average occupancy n — 0.8 and inverse temperature P — 8, 
16, and 30 eV"^ 

ened); i.e., the life time of the quasiparticles is very short. 

Since the temperature dep endence of the spectrum was 
already discussed in Ref. '51, we restrict ourselves here to 
the optical conductivity in Fig. 1121 The low energy opti- 
cal conductivity is suppressed as a result of the pseudo- 
gap in the spectral density. With decreasing tempera- 
ture, the optical spectrum is more suppressed since there 
are less thermal excitations and the two peak lattice dis- 
tortion becomes more pronounced (Fig. [7]) as does the 
pseudo gap^. Hence, the resistivity is enhanced at low 
temperature, giving rise to the insulating behavior in the 
paramagnetic phase, see Ref. [51| for the T-dependence of 
the resistivity. 

The two peaks of the optical spectrum stem from the 
two peaks in the electronic spectrum above ther Fermi 
energy. An excitation of an electron into the unoccupied, 
undistorted states just above the Fermi energy gives rise 
to the midgap states around 1.5 eV. Transitions to the 
distorted states far above the Fermi energy with two Cg 
electrons are at the origin of the peak around 6 eV. 



B. Ferromagentic phase transition 

The low temperature ferromagnetic phase is a bad 
metal which due to the strong spin and phonon scatter- 
ing has quasiparticles with very short life times and no 
true Drude peak^. Here, we will concentrate on the fer- 
romagnetic phase transition which can be either driven 
by applying a magnetic field or by decreasing the temper- 
ature. Fig. [T3l shows the doping dependence of the Curie 
temperature for the Jahn- Teller coupling g — 0.10 eV'^/^. 
Compared to the experimental results, the theoretical 
predictions of Tc have a similar shape with a maximum 
at intermediate doping but overall the values of Tc are 
about 2-3 times larger in magnitude. Close to the un- 
doped parent compound {x — 0, n — I), the Curie tem- 
perature is suppressed since the Coulob interaction hin- 
ders the movement of the electrons so that the double ex- 
change is no longer effective; instead antiferromagnetism 
prevails'*'^ . 

The magnitude of the Cure temperature depends 
strongly on the strength of the Jahn- Teller coupling and 
the hopping integral of the Cg electrons. Fig. [T3] plots its 
coupling dependence at n = 0.8, which are similar to the 
previous result^iS^. The Curie temperature decreases 
rapidly with increasing coupling. For strong Jahn- Teller 
coupling, a slight increase of g from 0.10 to 0.12 (or the 
dimensionless coupling A from 2.0 to 2.2) reduces Tc by 
a factor of 2 so that g = 0.12 eV^/^ (or A = 2.2) gives the 
correct experimental value. 

The Curie temperature may also be affected by many 
factors which are not included in our low-energy Hamil- 
tonian. These include the quantum fluctuation of Mn t2g 
spins and the antiferromagnctic supcrexchangc coupling 
between Mn t2g spins. Both tend to weaken the fer- 
romagnetic order and suppress the Curie temperature. 
The antiferromagnctic superexchange coupling has been 
estimated to be the order of 200 Ki^iS^ which, if taken 
into account, would greatly reduce the theoretical value 
of the Curie temperature. Also the mean field character 
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FIG. 14: Curie temperature as a function of the dimensionless 
coupling X = g/Q,y/to at n = 0.8. 



of the DMFT approximation tends to overestimate the 
Curie temperature. 

In the Hterature, the Curie temperature has been cal- 
culated with different methods such as DMFT^i^iia, 
conventional mean-field theory^, QMC simulations^, 
and the many-body CPA^**. Some of the results seem 
to be in better agreement with experiments. However, 
we should note that these results are all based on the de- 
tails of the models and approaches and are very sensitive 
to the values of the parameters which, unfortunately, are 
not always reliable and, as a matter of fact, vary con- 
siderably in the literature. A complete analysis of the 
problem is still required. 
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FIG. 15: Probability distribution of the breathing (Qi) and 
Jahn-Teller (Q) modes for the model Hamiltonian (Tj). The 
parameters are n = 0.8 and /3 = 16eV~^, U — 5eV, U' = 
3.5 eV, J = 0.75 eV, g = 0.08eV^/^ and = 0.07 eV for all 
three modes. 
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V. BREATHING MODE 

In this section we include a third electron-phonon cou- 
pling, the breathing mode which couples to the electron 
density, see section |IT1 Fig. [T51 shows the probability dis- 
tribution of both the breathing and Jahn-Teller phonon 
fields at /3= 16eV-\ n = 0.8, g = 0.08 eV^/^, [/ = 5eV, 
J = 0.75 eV and U' — 3.5 eV, and a phonon frequency 
il ~ 0.07 eV for all three modes. Except for the Jahn- 
Teller distortion Q being positive by definition, the distri- 
bution function of the breathing phonon is similar to that 
of the Jahn-Teller distortion and also has two peaks lo- 
cated at about and . These two peaks are related 
to the polaron states and the midgap states discussed be- 
fore. The only difference is the way in which the phonon 
modes are coupled to the Cg electrons: The breathing 
mode is coupled to the electron density, while the Jahn- 
Teller modes are coupled to the orbital polarization (with 
respect to a certain basis). Due to the Hund's couphng 
and the strong Coulomb interaction, double occupations 
are forbidden. Hence, the difference is not reflected in 
the distribution function. 

Since the breathing mode only couples to the electron 
density, it lowers the localization energy of the polaron 



-4 -2 2 4 6 8 10 
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FIG. 16: Electron spectral density for U' = 3.5 eV, n = 0.8 
and (3 — 16eV~^. The coupling constant is taken as g — 
0.08 eV^^^ for calculations with all three phonon modes and 
g = 0.08, 0.10 eV^/^ for calculations with only the Jahn-Teller 
modes. The phonon frequencies are taken as Q — 0.07 eV for 
all three modes. 



states but leaves the midgap states unchanged. The Cg 
electrons are thus more localized due to the inclusion of 
the breathing mode and the system becomes more insu- 
lating. 

Fig . [TBI compares the spectral densities at /3 = 16eV~^ 
and n — 0.8 with and without the breathing phonon. 
As expected, we see the density of states at the Fermi 
energy to be strongly suppressed by the existence of the 
breathing mode, demonstrating how the breathing mode 
supports the tendencies towards localization of the Cg 
electrons. 
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ln=1.0 
•n=0.825 



U(eV) 

FIG. 17: Two sets of parameters which fit the exper- 
imental energy gap in LaMnO^ and the resistivity in 
Lao.s^bSro.nsMnO^. Their intersection gives a single set 
(U,g ) which is expected to describe quantitatively both doped 
and undoped manganites. The insets show the DMFT (QMC) 
results as a function of the Coulomb interaction U and the 
Jahn-Teller coupling g for: (a) the energy gap at n = 1 
and (b) the resistivity at n = 0.825. Both are calculated 
at (3 = 30 eV"'^. The experimental results are indicated by 
the dotted line with the 1 eV energy gap in LaMnO^ and the 
resistivity of about 0.035 flcm for Lao,s25Sro.i75Mn03^. 



A. Determing U and g 

Let us now try to estimate from our results including 
the breathing mode the two parameters with the hitherto 
largest uncertaity, i.e., the on-site Coulomb interaction U 
and the electron-phonon coupling g. All other parame- 
ters are fixed as discussed in Section |lll To this end, we 
fit the two parameters to two experimental values, i.e., 
the experimental gap of undoped LaMnOs which is 1 
eV and the resitivity for x = 17.5% Sr doping which is 
0.03517cm. The calculations are done for (3 = 30 eV"^. 

The inset (a) of Fig. [T7] plots the energy gap for differ- 
ent Coulomb interactions U and Jahn-Teller couplings g 
calculated by DMFT (QMC) for the model Hamihonian 
As expected, the energy gap depends strongly on 
both parameters. 

The inset (b) of Fig. [l7] shows the DMFT (QMC) re- 
sults of the resistivity for different U and g at n — 0.825 
{x = 0.175) and /? = 30eV~^. As has been analyzed be- 
fore, the Coulomb interaction also affects the resistivity 
of the system since it enhances the electron localization. 
This is now proved explicitly in the figure. The resistiv- 
ity increases with increasing U, albeit it depends more 
sensitive on g than on U . 

By comparing with the experimental data (dotted lines 
in the two insets), the two experiments provide for two 
different sets of U and g parameters which we plot in 
the main panel of Fig. [T71 Only the intersection at 
U = 3.3 eV and g = 0.077 eV'^/^ agrees with both exper- 
iments. These parameters determined by us agree with 



w 


U 


2J|S| 


J 


Q. 


g 


3.6 eV 


3.3 eV 


2.7 eV 


0.75 eV 


0.07 eV 


0.077 eV^/^ 



TABLE I: Parameters estimated for doped and undopd man- 
ganites. W = 6to: bandwidth; U: intra-orbital Coulomb 
interaction; J: Cg-eg Hund's exchange; J: eg-t2g Hund's cou- 
pling; Q: phonon frequency; g: Jahn-Teller coupling. The 
bandwidth is obtained from the LDA calculations for the cu- 
bic structurelS, the Hund's coupling is calculated by the con- 
strained LDA for the ferromagnetic phase, and the phonon 
frequency is estimated from the Raman spectroscopy^^:. Only 
the Coulomb interaction U and the Jahn-Teller coupling g are 
estimated from the DMFT (QMC) calculations for the model 
Hamiltonian ^ by fitting the experimental data of the energy 
gap in LaMnOs and the resistivity in Lao.825Sro.i75Mn03, see 
Fig. [13 
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FIG. 18: Optical conductivity at doping x — 0, 0.1, 0.175 
in the paramagnetic (PM) phase and at x — 0.3 in the fer- 
romagnetic (FM) phase. The parameters are taken as is in 
Table\^and the inverse temperature is (3 = 30eV~^. Tiie in- 
set shows the effective carrier concentration calculated from 
the integrated optical conductivity. 



the crude estimates of Section |TT1 We hence expect this 
paremeter set, summarized in Table U to be the proper 
set for the model Hamiltonian ^ , describing both doped 
and undoped manganites. 



B. Optical conductivity 

The optical conductivity with the breathing mode in- 
cluded for the determined set of paramaters is presented 
in Fig. [181 The inset is the effective carrier concentration 
Neff{uj) calculated from the integral of the optical con- 
ductivity which can be compared directly to the experi- 
mental results of Ref. H^. We see a continuous crossover 
from the undoped to the doped system and a good agree- 
ment in the general behavior of the optical conductivity 
with experiment. However, the experimental magnitude 
of the optical peaks are 2-3 times larger than the the- 
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oretical one--. This gives rise to an unexpected large 
experimental Nf.ff{ijj) even below 4 eV, in contrast to 
the theoretical predictions presented in the inset of Fig. 
[T51 This suggests that, besides the tg states also oxy- 
gen p and Mn <2g states which we did not account for in 
the theoretical calculation contribute to the experimen- 
tal optical conductivity. Another possible explanation is 
the insufficiency of the group velocity as a substitute for 
the more apropriate dipole matrix elements. 

VI. CONCLUSION 

We have used DMFT to study the physics of mangan- 
ites by hands of a realistic microscopic model which takes 
into account both the electron-electron and electron- 
phonon interactions, together with the Hund's rule cou- 
pling between the Cg conduction electrons and the t2g 
spins. In the undoped system, the model produces simi- 
lar results as the previous LDA-f-DMFT calculations and, 
most surprisingly, it also predicts the correct structural 
transition temperature from dynamic to static Jahn- 
Teller distortion. In the doped phase, the eg electrons 
are trapped by the large lattice distortion and form a 
lattice polarons. This process is strongly supported by 
the Coulomb interaction. Our results provide for an ex- 
planation of the insulating-like paramagnetic state over a 



wide range of doping. The CMR is a result of a transition 
towards a ferromagnetic (bad) metallic state at a Curie 
temperature which shifts by applying an external ma- 
netic field. The dynamical properties are determined by 
the polaron states and midgap states which stem from the 
undistorted unoccupied states slightly above the Fermi 
eenrgy. The combination of both gives rise to the pseudo- 
gap behavior observed in doped manganites. The inclu- 
sion of the breathing mode further favors the tendencies 
towards polaron formation. Our results show that the 
realistic microscopic model can be applied to both doped 
and undoped manganites and can therefore be taken as 
the starting point towards a complete understanding of 
the physics of manganites in which the cooperative effect 
of the Jahn- Teller phonons and the quantum fluctuation 
and superexchange coupling of the t2g spins must also 
be taken into account for quantitative comparison with 
experiment. For future studies we provide for a realistic 
set of model parameters in Table HI 
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